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Abstract 

The chiral properties of the pseudoscalar mesons are studied numerically 
on a quenched 20^ lattice with the overlap fermion. We elucidate the role of 
the zero modes in the meson propagators, particularly that of the pseudoscalar 
meson. The non-perturbative renormalization constant Za is determined from 
the axial Ward identity and is found to be almost independent of the quark 
mass for the range of quark masses we study; this implies that the O(a^) error 
is small. The pion decay constant, /tt, is calculated from which we determine 
the lattice spacing to be 0.148 fm. We look for quenched chiral log in the 
pseudoscalar decay constants and the pseudoscalar masses and we find clear 
evidence for its presence. The chiral log parameter 6 is determined to be in the 
range 0.15 - 0.4 which is consistent with that predicted from quenched chiral 
perturbation theory. 



PACS numbers: 11.15.Ha, 12.38.Gc, ll.30.Rd 



1 Introduction 



One of the main goals of lattice QCD is to understand from first principles low- 
energy phenomenology as a consequence of chiral symmetry. Recent advances in the 
formulation of chiral fermions on the lattice hold great promise for studying chiral 
symmetry of QCD at finite lattice spacing 

Neuberger's overlap fermion 0], derived from the overlap formalism [Q], is such 
a chiral fermion on the lattice and has been implemented numerically to study the 
chiral condensate [|, H, 0, quark mass H, 0, renormalization constants H, 0, 0, 
and short- distance current correlators and to check chiral symmetry [Til and 
scaling |^. However these studies are limited to small volumes due to the large 
numerical cost associated with approximating the matrix sign function. In this paper, 
we shall study physical observables, such as the pseudoscalar meson masses and pion 
decay constants, close to the physical u, d quark mass. As such, we need to work on a 
lattice which is at least 3 times larger than the Compton wavelength of the pion with 
the smallest mass in order to alleviate finite volume effects. We work on a 20^ lattice 
with a = 0.148 fm as determined from the pion decay constant /^r. This gives the 
lattice size La = 3.0 fm and the smallest pion mass is ~ 280 MeV. Thus, the lattice 
size is ~ 4 times the Compton wavelength of the lowest-mass pion. 

This paper is organized as follows. We will give the numerical details of the calcu- 
lation in Sec. 2. In Sec. 3, we shall discuss the effect of the zero modes in the meson 
propagators. In view of the fact that the scalar condensate receives a contribution 
from the zero modes (which goes away in the infinite volume limit) through the gener- 
alized Gell-Mann-Oakes-Renner (GOR) relation, the pseudoscalar correlator should 
also be contaminated by the zero modes. We have observed the effect of the zero 
modes in the pseudoscalar propagator at small quark mass. After clarifying the zero 
mode issue, we proceed to calculate the non-perturbatively determined renormaliza- 
tion constant Za from the axial Ward identity and the pion decay constant We 
find that /tt is free of the quenched chiral log singularity and has a small error in the 
chiral limit. Thus, it is a good quantity to set the lattice scale. We present the results 
in Sec. 4. In Sec. 5, we explain our effort in searching for the predicted quenched chi- 
ral logs. We see the chiral logs in the pseudoscalar masses, the pseudoscalar matrix 
element fp, and the fp/fn ratio at very small quark masses. A summary is given in 
Sec. 6. 
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2 Numerical Details 



For Neuberger's overlap fermion p[, we adopt the following form for the massive Dirac 
operator |, |g, |T3[ 

D(mo) = (1 - '^)pD{p) + moa, (1) 

where 

D{p) = l + ^,e{H), (2) 

so that 

^(-o) = P+^ + (p-^)75e(n (3) 

where e(iJ) = iJ/ V if ^ is the matrix sign function and H is taken to be the hermitian 
Wilson-Dirac operator, i.e. H = Here Dy^ is the usual Wilson fermion operator, 

except with a negative mass parameter —p = 1 /2k, — 4 in which < n < 0.25. We 
take K = 0.19 in our calculation which corresponds to p = 1.368. The massive overlap 
action is so defined so that the tree-level renormalization of mass and wavefunction 
is unity. The bare mass parameter [|, mo, is proportional to the quark mass without 
an additive constant which we have verified numerically in a previous study [Q. 

We adopt the optimal rational approximation [|^, |^ to approximate the ma- 
trix sign function. The inversion of the quark matrix involves nested do loops in this 
approximation. It is found that it is cost effective to project out a relatively few eigen- 
modes with very small eigenvalues in the operator in order to reduce the condition 
number and speed up the convergence in the inner do loop |0, §. At the same time, 
this improves chiral symmetry relations such as the Gell-Mann-Oakes-Renner rela- 



tion p. However, it is shown [15| that the density of these small eigenmodes grows 
as with a being the lattice spacing. As a result, it is very costly and impractical 
to work on large volumes with the lattice spacings currently used. There are simply 
too many small eigenmodes to be projected out. 

For this reason, we explore other options to clear this hurdle. We have tested the 



tree-level tadpole-improved Liischer-Weisz gauge action [|I6[ and find that the density 
of these small eigenvalue modes is decreased to a point where it becomes feasible 
to go to large volumes with a lattice size 4 times the Compton wavelength of the 
lightest pion. We further find that the anisotropic action |]T^ requires projection of 
more small eigenvalues in in order to achieve the same convergence in the inner 
loop than does the isotropic one. Thus, we decide to use the isotropic action. We 
also find that using the clover action with either sign requires the projection of more 
small eigenvalue modes. Therefore we use the Wilson action for H in the Neuberger 
operator. On a 20^ lattice with j3 = 7.60 tree-level tadpole-improved Liischer-Weisz 
gauge action, we project out 85 small eigenmodes. Beyond these eigenmodes, the 

*Notc that we used a different normalization in the action before in As a result the bare mass 
here is equal to p times the bare mass in 
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level density becomes large. As a result, the number of conjugate gradient steps is 
about 345 for the inner loop and about 300 for the outer loop. While the number for 
the inner do loop seems to be fairly independent of the lattice volume, the number 
for the outer do loop is about a factor of 2 larger than those for the Wilson gauge 
action on small volumes P). 

Since the conjugate gradient algorithm accommodates multiple masses with a 
minimum overhead, we calculated 16 quark masses ranging from rriQa = 0.01505 to 
■rriQa = 0.2736 which are listed in Table 1 together with the number of configurations 
for each mass. 

From the string tension with ^/o^ = 440 MeV, we find that a = 0.13 fm. However, 
as we shall see later in Sec. 4, the scale determined from is 0.148(2) fm which 
makes the physical length of the lattice to be 3.0 fm. The smallest pion mass turns 
out to be ~ 280 MeV so that the size of the lattice is ~ 4.1 times of the Compton 
wavelength of the lowest mass pion and more than 4 times for the heavier ones. 

We adopt the periodic boundary condition for the spatial dimensions and the fixed 
boundary condition in the time direction so that we can have effectively a longer range 
of time separation between the source and sink to examine the meson propagators 
with small quark masses which lead to longer correlation lengths. The source of the 
meson interpolation field is placed at the 3rd time slice and we consider the sink as far 
as the 16th time slice to mitigate the boundary effect. This gives us a time separation 
of 13. 

We have varied statistics for different quark masses. Of the 16 cases, 7 have 25 
gauge configurations, one has 53, and the remaining 8 have 63 configurations. In order 
to carry out correlated fits to extrapolate observables to the physical pion mass, we 
construct the covariance matrix by embedding the one with smaller dimension, e.g. 
25 and 53 into the one with dimension 63 in a block diagonal form. For example, 
the covariance Cy for the one with 25 configurations is constructed so that Cij{i < 
25 J > 25) = 0, C{i > 25, j < 25) = and Cij{i,j > 25) = 6ij. 



Table 1: Quark mass rriQa and number of gauge configurations are listed. 
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0.05472 


0.06840 


0.08208 
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63 


53 


63 


25 
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Figure 1: {ipil)) as a function of the quark mass. We used 50 configurations of a 6^ x 12 
lattice with Wilson gauge action at /3 = 5.7. Here /ia = m^ajlp. 

3 Zero Mode Effects in Meson Propagators 

The quark zero mode is known to contribute to the vacuum scalar density (^/'V') on a 
finite volume. The latter can be written in the following form for small quark mass 
moFI 



nioV 



+ Co + cimo, 



(4) 



where Q is the topological charge which, according to the Atiya-Singer theorem, is 
the difference between the number of left-handed and right-handed zero modes (i.e. 
Q = ri- —n+) and has been shown to hold for overlap fermions or other local fermion 
actions which satisfy the Ginsparg- Wilson relation [|18], |19|. Since {\Q\) grows as VV, 
the zero mode contribution vanishes in the infinite volume limit while keeping rriQ 
fixed at a non-zero value. Thus the quark condensate which is the infinite volume 
and zero mass limit of (ipip) is represented by Cq in Eq. (^. However, on a fixed finite 
volume lattice, this zero mode contribution is divergent for small enough rriQ. This 
was first observed in the domain-wall formulation |^ and is also seen in the overlap 
fermion [^. Here we reproduce it in Fig. 1 which shows the divergent part of (ijjip) 
from the zero modes for a 6^ x 12 lattice with the Wilson gauge action at /3 = 5.7. 

We see from Fig. 1 that cq is non-zero in this range of the quark mass and upon 
extrapolation to the infinite volume before taking the chiral limit defines the quark 
condensate — S. However, if one keeps volume fixed and lets the quark mass approach 
zero, e.g. rriQa < 0.001, it is then found ^ ^ that cq becomes zero. It is known |^ 
that when the size of the lattice is much smaller than the pion Compton wavelength. 



^We shall address the quenched chiral log issue separately in Sec. |[ 
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2^* with 63 configurations ttloCl — 0.01915 




(t + to)/a 



Figure 2: Pion propagator for moct 
placed at to/ct = 3 



0.01915 for all 63 configurations. The source is 



i.e. L ^ I/jtItt, the constant term vanishes and (^/'^/') is proportional to mQT?V for 
small masses aside from the -^^^ term. Using finite size scaling, the chiral condensate 
S can be extracted [Q. 

While we have a reasonably good understanding of the role of zero modes in 
{ipip) , their role in the hadron propagators is only beginning to be investigated in the 
domain-wall formalism [^] and the overlap formalism |^ and its detailed influence 
on the hadron propagators is not fully understood. We shall investigate it in the 
pseudoscalar meson channel. There has been some concern about the behavior of the 
pion mass. It is not clear if it approaches zero in a finite volume in the chiral limit. In 
fact, one finds that the pion mass squared appears to approach a non-zero value on 
small lattices El, |23| . To understand the situation, we first examine the generalized 



Gell-Mann-Oakes-Renner (GOR) relation 



2(2mo)- 
V 



<fx(fy{'K\x)'K''{y)) = -2mo{ipip). 



(5) 



The left hand side of Eq. (13) has a contribution of 



2(IQI) 

V 



as does the right hand 



side. Assuming the remainder of the pseudoscalar susceptibility is dominated by the 
pion, it is approximately — /^m^. Comparing with Eq. (|), we get 



f^ml = 2moCo + 2mlci. 



(6) 
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20"* with 19 Q = configurations moa = 0.01915 

10 F 1 1 1 r- 



5 -r 



(t + to)/a 



Figure 3: Pion propagator for mga = 0.01914 for 19 configurations with Q = 0. 

From this we see that approaches zero as mo for those volumes where cq is 
non-zero. On the other hand, when the mo is so small that cq becomes zero for certain 
fixed volume, the r.h.s. of Eq. (|]) is survived by the next leading term with ci which 
is proportional to mg. As a result, one expects that m^r to be linearly proportional 
to mo- Here we have ignored the complication due to the quenched chiral log which 
we will address in Sec. 5. 

Next, we turn to the zero-momentum pseudoscalar propagator / d^x{'K°'{x)'K°'{0)) 
which has two terms due to the zero modes as pointed out in the study of domain 
wall fermions 



J d^x{7r{x)7r{0)) = I d^x[ J2 



i,j=zeromodes ^0 

^ 2 y tr{iji{x)Ux))tr{ijl{0)MO)) ^ , l(Ok(0)k)Pe— 
i=o,A>o mo(A2 + mg) 2m^ 

The first term is purely the zero-mode contribution. The second term is the cross 
term between the zero modes and the non-zero modes. We have used the property 
that the non-zero modes come in pairs which are related by 75, i.e. j^i/jx = ip-x- 
Upon integrating the propagator with respect to time, we find 



J d^x{TT{x)7T{0)) 



mo mo(A2 + mo) 2ml 

(8) 
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20"* with 44 Q 7^ configurations moa = 0.01915 

10 F 1 1 1 r- 
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Figure 4: The same as in Fig. 3 for 44 configurations with Q 0. 

Comparing with the generahzed GOR relation in Eq. (^), we see the first term corre- 
sponds to term in Eq. (^) and the second term vanishes due to the orthogonahty 
between the zero modes and the non-zero modes. In either case, we expect that the 
number of zero modes grows with and the eigenf unction ■4>{0) oc l/y/V. As a 
result, both zero- mode terms in Eqs. (|^) and (|[) decrease with volume like l/y/V 
and are finite volume artifacts. 

The most straight-forward way of isolating the zero mode contribution is to cal- 



culate their eigenvectors and project them out as is done in Ref. [^]. However, this 
is very costly, as costly as calculating the quark propagator itself. 

We show the pseudoscalar propagator for a light quark mass (moa = 0.01915) in 
Fig. 2 and we see that there appears a kink att/a~8 — 9 (The source is placed at 
to/a = 3 so that it appears at {t + tQ)/a 11 — 12.) To further explore its origin, we 
separate the 63 configurations into 19 with trivial topology (i.e. Q = 0) and 44 with 
non-trivial topology (i.e. Q 0) and plot the respective pseudoscalar propagators in 
Figs. 3 and 4. 

We see that the propagator on the Q = configurations (Fig. 3) has a single 
exponential all the way from t/a = 5 to 14. Upon fitting a single exponential in this 
range, we find mpa = 0.153(12). On the other hand, the propagator of the Q 
configurations (Fig. 4) still has a pronounced kink at t/a ~ 8. When we fit it in the 
range t/a = 8 — 14, the mass is mpa = 0.146(32) which is quite consistent with that 
from the Q = configurations. On the other hand, when we fit the time range t/a 
from 4 to 8, mpa = 0.245(58) which is quite a bit higher than that of the Q = 
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20^ lattice with Overlap Fermions t/a=5-Q 




O.OB 0.1 0.15 0.2 



rriQa 

Figure 5: mj^a^ vs rriQa with a fit of mpa in the window oft/a = 5 — 9. The dotted 
fine is a fit hnear and quadratic in mga. The sohd fine is a fit including the chiral log. 



configurations. We take this as the evidence that the zero-mode contributions, the 
combined direct and cross terms, fall off faster in t/a than the pseudoscalar mass. 
We can thus use the time separation as the filter to obtain the masses and decay 
constants of the physical pseudoscalar mesons. Now we can understand why in the 
previous studies of nip, the does not approach zero with the quark mass 
In Ref. p^], the lattice size is 6"^ x 12 with /3 = 5.7. Since it uses the fixed boundary 



condition in the time direction, the maximum usable time separation is about 7. This 
translates into time separation of ~ 9 on our lattice. In Ref. [|25[, the lattice size is 
12^ X 24 with /5 = 5.9. In this study, the authors use the periodic boundary condition 
in time, time slices up to 12 are fitted which corresponds to ~ 10 on our lattice. In 
either case, the fitted time range is expected to be contaminated by the zero mode 
contribution and result in a higher mass for a given quark mass. To verify this, we 
fit our data in the time range t/a = 5 — 9 and plot the resulting m^pO? in Fig. 5. 
We see that they indeed don't approach zero with a fit linear and quadratic in rriQa 
(dashed line), similar to those shown in Refs. |24, 2^. One may attempt to interpret 
the data to include a chiral log and force the pion mass to go to zero in the chiral 
limit. This can be misleading. We shall defer the discussion of the complication due 
to the quenched chiral log in Sec. 

In principle, one can overcome the zero mode problem by fitting the zero mo- 
mentum pion propagator Gpp{'p = 0,t) in a time range beyond the localized zero 
mode contribution, such as beyond {t + to)/a > 12 in Fig. 2. Unfortunately, our 
data at larger time slices are tainted by the refiection from the fixed boundary at 
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{t + t{))/a = 20 when the quark mass is small so that the fit in the limited time win- 
dow to avoid both the zero mode and the boundary reflection is unsettling in this case. 
For this reason, we examine the propagator Ga4,p{p = 0,t) instead. Since the zero 
modes on one gauge configuration have the same chirality, the pure zero mode contri- 
bution (i.e. the direct zero mode contribution which corresponds to the first term on 
the right- hand- side of Eq. (J^)) vanishes. The cross-term between the zero modes and 
the non-zero modes does not vanish, but is expected to be small due to cancellations. 
We plot in Fig. 6 GA^pip = 0,t) for the same quark mass niQa = 0.01915 and we do 
not see a kink in the range 5 < {t + to)/a < 15 a.s is in Gpp{p = 0,t) in Fig. 2. Thus, 
we shall calculate the pseudoscalar masses from GA^pip = 0, t). The results of mpO^ 
so obtained are plotted in Fig. 7. 

Comparing with Fig. 5, we note that m^pO^ in Fig. 7 which are free from the zero 
mode contamination are lower than the corresponding ones in Fig. 5 for small quark 
masses. How m?pa^ approaches zero in the chiral limit is complicated by the presence 
of the quenched chiral log. We shall postpone this discussion until Sec. |^, except to 
mention that the solid line in Fig. 7 is the fit with the chiral log and the dashed line 
is the fit with linear and quadratic terms in niQa. 



63 configurations with mo = 0.01915 



^ 0.001 

e3 



(t -+- tQ)/cL 



Figure 6: Propagator Ga4,p{p = 0,t) for moa = 0.01915 for all 63 configurations. 
The zero mode contributions to mesons can be written as 

d^x{M{x)M{Q))\zeromodes= (9) 
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0.6 




rriQa 

Figure 7: m^pO^ vs rriQa calculated from Ga4p{p = 0,t) The linear plus quadratic fit 
(dotted line) and the chiral log fit (solid line) from Eq. (|32D with A^^ = 0.8 GeV are 
discussed in Sec. ^ 

tr{^lJUxh,T^|J,{x))tr{^pl{0)T^,^}Jj{0)) 



i,j=zero modes 



ml 



-2 E 

i=0,A>0 



mo(A2 + ml) 



(10) 



where F and F are the gamma matrices for the corresponding meson interpolation 
fields. For the pseudoscalar meson as we discussed above, F = — f = 75. For the 
scalar meson (the connected insertion part), F = F = 1. Since the zero modes are 
the eigenstates of 75, i.e. 75'?/'j=o = ±'?/'j=05 the zero mode contribution in the scalar 
propagator has a negative sign from that in the pseudoscalar propagator. Thus, it is 
suggested [^, H, 0] to consider / (i'^x[(7r(x)7r(0)) + ((T(x)cr(O))] where the zero-mode 
contributions cancel and at large t/a it should be dominated by the pseudoscalar. 
However, there is a problem. As will be shown in a separate publication |2^, the 
large time part of the isovector-scalar propagator for the quark mass range that we are 
concerned turns out to be negative. It is pointed out in a study with pole-shifting in 
the Wilson action that it is dominated by the would-be 77' and tt intermediate state 



which is negative due to quenching. In the intermediate time region, there is still the 
contribution from rj' and vr intermediate state besides the Qq. As a consequence, the 
addition of the pseudoscalar and scalar meson is not a viable solution to obtaining the 
pseudoscalar mass. It is also suggested that the axial-vector interpolation field with 
F = F = 7475 does not have the direct term contribution from the zero modes since all 
the zero modes have the same chirality in a given gauge configuration. However, the 
second term in Eq. (|^) may still have a contribution. Since it is a cross term, it might 
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be small due to cancellations. Unfortunately, our data on the (AiA^) correlator are 
much noisier than the (tttt) correlator, since (0|y44|7r(0)) = v^^ti^/tt which diminishes 
when the pion mass is small. We cannot conclude anything from them. 

It appears that, short of projecting out the zero modes from the meson propaga- 
tors, fitting the pseudoscalar mass from the propagator GAiPip = 0, t) is probably the 
only practical way of getting reasonably reliable and accurate pseudoscalar masses 
for our lattice with fixed boundary condition. 



4 Za and Pion Decay Constant /^r 

It has been pointed out in our earlier work that the renormalization constant Za 
for the axial current = ip^i'y^'y^^l — D/2p)^)ip can be obtained directly through 
the axial Ward identity 

ZaO^A^ = 2Z^moZpP, (11) 

where P = — D/2p)^)ip is the pseudoscalar density. Since Z^ = Z~g" and 

Zs = Zp due to the fact that the scalar density ipll — D/2p)^ip and the pseudoscalar 
density P are in the same chiral multiplet, Z^ and Zp cancel in Eq. (|TTp and one 
can determine Za to O(a^) non-perturbatively from the axial Ward identity using 
the bare mass rriQ and bare operator P. To obtain Za, we shall consider the on-shell 
matrix elements between the vacuum and the zero-momentum pion state for the axial 
Ward identity 

ZA{0\d^A^\7iip = 0)) = 2mo(0|P|7r(p = 0)), (12) 
where the matrix elements can be obtained from the zero-momentum correlators 

Ga,AMP=0^t) = {J2d,A,{x)PiO)) 

X 

Gpp{p=0,t) = (5:P(x)P(0)). (13) 

X 

The non-perturbative Za is then 

Za = hm ^rnoGpp{p = 0,t) 



There has been extensive study 30, Q of the 0(a) improvement of the 



Wilson action and composite operators in relation to the chiral symmetry and axial 
Ward identity. It has been shown that in the improved mass-independent renor- 



malization scheme, the renormalized improved axial current and pseudoscalar density 
have the following form from an 0(a) improved action 



A^ = ZA{l+bAmga){A^ + CAad^P}, 

P^ = Zp{l + bpmga)P, (15) 
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20^ Lattice with Overlap Fermions 




rriQa 

Figure 8: vs quark mass moO. The results from Eq. (0) are indicated by o and 
those from Eq. ([T7|) are indicated by •. 



where rriq = mQ — nic is the subtracted quark mass and ca-, and hp are improvement 
coefficients. The renormahzation constants Za and Zp are functions of the modified 
couphng = (7q(1 + hgrnqo). Now that the overlap fermion is 0(a) improved [p3| , 



it satisfies the Ginsparg- Wilson relation and the quark mass is not additively 



renormalized. As a result, hA = ca = hp = 0. The renormahzation constant Za 
determined from Eq. (0) has only 0{o?) error as a consequence. Za thus computed 
for all the quark masses are plotted in Fig. 8 as indicated by the open circles. We see 
that there is still an appreciable 0{a?) error. Given that the time derivative itself in 
GdiAiPip = 0,t) invokes an 0{a?) error, it would be better to adopt a definition for 
Za which is devoid of this superfiuous O(a^) error. This can be achieved by noticing 
that, at large t where the pion state dominates the propagator Gq^a4,p{p = 0,t), one 
can effectively make the substitution 

GdiAip{'P = 0,t)^"- m^GAipip = O.^)- (16) 



Consequently, Eq. (|14|) becomes 



2moGpp{p = 0,t) 

^A = Jim — — — -. (17) 

m^GAiP[p = 0,t) 



t- 



We plot the results of Za from Eq. {\V7\} in Fig. 8. Save for the last two points at 
the smallest masses (moa = 0.01505 and 0.01642), the errors are small. We should 
point out that it is conspicuously flat as a function of mga indicating that the O(a^) 
error from the action and the operators is small. We flt them in the form Za + hmla? 
and found that Za = 1.589(4), 6 = 0.175(73) with x^/DF = 0.30. We see that Z 
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Table 2: Renormalization constant Za from Eq. (|T^ (column 2), and Eq. (|T^ 
(column 3). 





Za (Eq. (0)) 


Za (Eq. (0)) 


0.2736 


1.490(2) 


1.605(4) 


0.2462 


1.499(3) 


1.594(6) 


0.2189 


1.511(2) 


1.597(4) 


0.1915 


1.519(3) 


1.589(8) 


0.1642 


1.532(2) 


1.596(5) 


0.1368 


1.539(4) 


1.585(11) 


0.1094 


1.551(2) 


1.592(8) 


0.09576 


1.552(4) 


1.587(12) 


0.08208 


1.560(2) 


1.593(8) 


0.06840 


1.561(5) 


1.594(13) 


0.05472 


1.569(3) 


1.586(12) 


0.04104 


1.572(5) 


1.595(15) 


0.02736 


1.576(8) 


1.604(17) 


0.01915 


1.583(17) 


1.585(20) 


0.01642 


1.585(72) 


1.618(60) 


0.01505 


1.596(105) 


1.598(114) 



is determined to the precision of a fraction of 1%. We find it to be larger than the 
perturbative calculation |12| which gives .^^(At = 1/a) = 1.213 for Wilson gauge action 
with (5 = 5.85 which has about the same string tension as our gauge configurations. 
Za from both Eqs. (|14D and ([TTD are tabulated in Table 2. 

For the pion decay constant, one can look at the ratio of the zero- momentum 
correlator 

GaAp = M = {Y^A,{x)PiO)), (18) 



and Gppij) = 0, t), 



Ua = hm ZAGAAP-O.t) ^^g^ 



t >oo 



■mT,a Gppip = 0, t) 



For our definition of the isovector axial and pseudoscalar currents, the experimental 
U is 92.4 MeV. 

Combining with Za from Eq. (|l^ , we obtain 



2moaJGpp{p = 0,t)GAAp = 0,t) 

Ua = hm " e"-*/l (20) 

t^oo y/m^aGQ^A4P{P = 0, t) 
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We first fit the pion masses from Ga4,p{'P = 0,t) and feed tliem into Eq. (|20|). 
Unlike tlie case witli in Eqs. (|l^) and ([TTD wliere tlie boundary effects on tlie 
common source for tlie interpolation field P in the numerator and denominator cancel, 
they don't cancel in /^r. To correct for this, we replace \jGpp{p = 0, t) with Gpp{p = 

0,t) e^-^l'^I^Gpp{p = 0, (iVt + l)a - to) e™-((^*+^)'^-*o) which should get rid of the 
boundary effect which affects the matrix element (0|P|7r) associated with the source 
in Gpp. The errors of the ratio are obtained with the jackknife method. We plot the 
result in Fig. 9 as a function of m^a as open circles. We also tabulate them in Table 
3 together with the pion mass computed from Ga^p{p = 0,t). 



20^ Lattice with Overlap Fermions 

0.12 I 1 1 1 , 1 




0.05 0.1 0.15 0.2 0.25 0.3 

rriQa 

Figure 9: Renormalized fj^a vs quark mass rriQa. 

As in the case of Za defined from Gd^A^pip = 0, t), /tt defined from Eq. (^) con- 
tains O(a^) error invoked by the time derivative. We shall again make the substitution 
in Eq. (|l6]) and arrive at 



2moaJGpp('p = 0,t) m^ae™''*/^ 
Ua = lim ^ "^"^l ; . (21) 



The results are listed in Table 3 and plotted in Fig. 9 as •. The difference from 
the o reflects the 0{a?) error due to the time derivative in Gs^AiPip = 0,t). With all 



14 



16 data points in a linear fit, we find 



f^a = 0.0691(11) + 0.109(56) niQa, (22) 

with x^/DF = 0.03. We note that it yields an error of only 1.6% in the chiral limit. 
Given that it is a physical observable from the fermion operator without the chiral 
log in a quenched theory and with such high precision, it is an ideal quantity to 



set the scale of the lattice. Comparing with the experimental value = 92.4 MeV, 
we determine the scale of our lattice to be a = 0.148(2) fm. This is about 14% 
higher than that determined from the string tension with ^/a = 440 MeV and 20% 
larger than that determined from tq. This is quite consistent with other quenched 
calculations. 

Our results on and its quark mass dependence are consistent with those calcu- 
lated with the domain-wall fermion |^ and the overlap fermion ^ on smaller lattices 
and higher pion masses. 



Table 3: Pion masses obtained from Ga4,p{p = 0, t) are listed together with f^^a from 
Eq. (gop (column 4), and Eq. (|21]) (column 5), and fpa^ from Eq. (g). 



iTiQa 






(Eq. (|0D) 


(Eq. (H) 


fpa' 


0.2736 


0.705(4) 


0.497(6) 


0.0925(10) 


0.0991(10) 


0.1273(13) 


0.2462 


0.660(7) 


0.436(9) 


0.0896(23) 


0.0952(23) 


0.1192(28) 


0.2189 


0.620(4) 


0.384(5) 


0.0883(12) 


0.0929(11) 


0.1155(14) 


0.1915 


0.573(8) 


0.328(9) 


0.0852(24) 


0.0893(31) 


0.1082(38) 


0.1642 


0.529(7) 


0.280(7) 


0.0839(14) 


0.0869(13) 


0.1048(16) 


0.1368 


0.477(15) 


0.228(14) 


0.0813(33) 


0.0837(33) 


0.0983(39) 


0.1094 


0.429(8) 


0.184(7) 


0.0794(13) 


0.0810(19) 


0.0963(22) 


0.09576 


0.399(17) 


0.159(14) 


0.0784(28) 


0.0794(37) 


0.0934(43) 


0.08208 


0.374(9) 


0.140(7) 


0.0770(13) 


0.0779(21) 


0.0939(25) 


0.06840 


0.343(19) 


0.118(13) 


0.0755(26) 


0.0762(40) 


0.0927(49) 


0.05472 


0.312(9) 


0.0973(56) 


0.0744(15) 


0.0753(18) 


0.0942(30) 


0.04104 


0.272(12) 


0.0740(65) 


0.0729(18) 


0.0733(25) 


0.0933(32) 


0.02736 


0.236(17) 


0.0557(80) 


0.0722(19) 


0.0728(29) 


0.1055(34) 


0.01915 


0.216(16) 


0.0467(69) 


0.0721(26) 


0.0709(27) 


0.1222(42) 


0.01642 


0.218(47) 


0.0475(205) 


0.0714(44) 


0.0708(58) 


0.145(12) 


0.01505 


0.209(53) 


0.0437(221) 


0.0734(50) 


0.0738(68) 


0.152(14) 
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5 Quenched Chiral Logs 



In the quenched approximation of QCD, one ignores the virtual quark loops. One 
of the consequences is that the flavor-singlet meson propagator (which we shall refer 
to as 77' even though the physical is not completely flavor-singlet and has octet 
mixture) has a double pole in the Veneziano model for the U(l) anomaly p7 |. 



As such, it does not move the mass of the would-be Goldstone boson to the large 77' 
mass. Consequently, this leads to infrared singular v( loops with the hair-pin type 
diagrams in the renormalization of hadron masses and certain matrix elements and 
thus alters their chiral behaviors from those of full QCD with dynamical fermions. 

The flrst study of the anomalous chiral behavior was done by Sharpe |^ and 
Bernard and Golterman in quenched chiral perturbation theory. They predicted 



the chiral-log pathologies in the pseudoscalar meson masses, ('?/'^/'), the JkI f-K ratio, 
etc. The flrst evidence of the chiral-log was observed by the CP-PACS Collabora- 
tion [PU| in the ratio of pseudoscalar meson masses with two unequal quark masses 



with the Wilson fermion. They obtain the chiral log parameter b = 0.8 — 1.2. A 
more extensive study [ED| which invokes the shifting of the real poles of the quark 



propagator to improve the otherwise poor chiral properties of the Wilson fermion 
near the chiral limit was carried out to examine the quenched chiral logs in the pseu- 
doscalar masses and the pseudoscalar decay constants, obtaining consistent results 
with 6 = 0.065 ± 0.013. The same quenched chiral log is also observed in the ratio of 
rn^/mq with the Kogut-Susskind action [|T| with 5 = 0.061 ± 0.03. All of them are 



small compared to that expected from the coupling of the would-be Goldstone bosons 
(or quark loops in the pseudoscalar channel) which is responsible for the rj' mass of the 
U(l) anomaly. Recently, the small nonzero eigenvalues of the Overlap Dirac operator 
were calculated in a deconflned phase |^ and it was found that the chiral condensate 



diverges at the inflnite volume limit indicating a quenched singularity consistent with 
the quenched chiral perturbation prediction (tptp) oc m~ . 

Given that the overlap fermion has the promise of exact chiral symmetry on the 
lattice, it is natural to look for these chiral singularities and check if the quenched 
chiral logs seen in the Wilson and Kogut-Susskind fermions can be verifled with the 
overlap action. The flrst attempt to extract the chiral log from the pion mass on 
several small volumes was inconclusive We shall examine them here on a much 
larger volume. 



The behavior of the quenched chiral logs can be seen from the sigma model |^| 
with ^7(3) X ^7(3) where the pseudoscalar fleld is represented by 

= g0o//giELi^<"^-//, (23) 

where Aa is the SU(3) flavor matrix and (pa are the octet Goldstone boson flelds. The 



U(l) part is described by r]' fleld with 0o = J'^/^fV'/ f where Nf is the number of 
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flavors which is 3 in our case. The effect of the chiral logs can be understood as the 
renormalization effect of integrating out the 77' |^0|. The resulting SU{3) x SU{3) will 
be represented by the renormalized U 

U = e-^'^o)/V'eELr>^'^'t>'^/f. (24) 

In the quenched approximation, the integral representing the 77' loop involves only 
the hair-pin diagram of two would-be singlet Goldstone boson (we shall refer to it as 
rj and it has the same mass as vr) propagators 



Therefore, the infrared singular part of U in Eq. ( p4D can be represented by 6 
where 



'ml 



From the Witten-Veneziano model of the 77' mass, rfiQ ~ 870MeV. This gives an 
estimate of 5 = 0.183. 

To see the effects on various physical quantities, one can first look at pseudoscalar 
density and axial current operators |?D|. In the sigma model, 

ijji'y^il) oc U — U'^ 
^27^75^' « i[U~^d^U -{d^U~^)U]. (28) 

With U given in Eq. (|26|), one arrives at 

1 



m 



p 

U = (0|V'^7M75V^k(p = 0))/V2m^ = ^ (29) 

where fp is the pseudoscalar decay constant and /^r is the axial decay constant and 
fp and are constants for small quark masses. Thus, one expects that that fp is 
singular as the quark masses approaches zero in the quenched approximation; whereas 
remains a constant. From the axial Ward identity in Eq. ( [TT| ) and Eq. (^), one 
expects 

2 fp I ^ \& 
ml 

where mq is the quark mass and therefore 

,2 ^TTJ 



ml oc mq— oc mq{—Y, (30) 



m; oc ml+' (31) 
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which is the behavior predicted in quenched xPT [^. 

Given m^a fitted from Ga4p{p = 0,t) and hsted in Table 3, we first fit them in 
the form in Eq. (|31|) for a range of small quark mass points. We list the fitted 6 and 
x'^/DF in Table 4. We find that S ranges from 0.145(66) to 0.24(12). Beyond this 
range of 8 to 11 smallest quark masses, the errors are greater than half of the fitted 
6 value. 

Table 4: Quenched chiral log parameter 6 and x^/DF as fitted from irn^o? in Eq. 



^ of smallest rriQa 


5 


XVDF 


8 


0.24 (12) 


0.24 


9 


0.224(91) 


0.22 


10 


0.163(75) 


0.28 


11 


0.145(66) 


0.28 



Since the form in Eq. (pT]) is limited to small mga, we shall also fit m^a^ with the 



form |8|, |9 



mla^ = Amoa{l - 6[\n{Amoa/Ala^) + 1]} + Em^a^ (32) 

which allows a fit to cover the whole range of 16 quark masses. The best fits which 
give stable values of A and B and with errors less than half of the fitted values of 6 
for a range of A^^ = 0.6 GeV to 1.4 GeV are listed in Table 5. 

Table 5: Quenched chiral log parameter 5 and / DF as fitted from in Eq. (R^). 



Ax (GeV) 


A 


B 


5 


XVDF 


0.6 


1.72(7) 


3.0(8) 


0.23(7) 


0.18 


0.8 


1.42(11) 


3.0(8) 


0.28(11) 


0.18 


1.0 


1.17(24) 


3.0(8) 


0.34(17) 


0.18 



They are consistent with those fitted in the exponential form, albeit with values 
of 5 between 0.2 and 0.3 which are somewhat higher than those from the exponential 
form. 

We plot the fit with A^ = 0.8 GeV as a solid line in Fig. 7. To check if there is 
indeed a bona fide quenched chiral log, we fit rn^a? alternatively with powers in niQa 



18 



20^ Lattice with Overlap Fermions 

0-3 I 1 , 1 1 , 



0.25 - 



0.2 
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Figure 10: Renormalized fpo^ vs quark mass mga. The solid line is a fit excluding 
the two smallest quark masses with A^^. = 0.8 GeV in Eq. (|35D . 

up to cubic terms and found 

mla^ = 1.65(5) moa + 0.54(21) m^a^, x'^/DF = 0.90, 

mlc? = 1.88(10) moa- 2.3(11) m2a2 + 7.8(30) m|^a^ x^/DF = 0A6. (33) 

The fit including the cubic term in mga is shown in Fig. 7 as the dotted line. We 
see that in either case, the / DF is larger than that with the quenched chiral log in 
Table 5 which is 0.18. From this, we conclude that the quenched chiral is seen in our 
data for m^a^. 

We should point out that if we were to use the m?pO? data from the shorter time 
range such as the one with t/a = 5 — 9as plotted in Fig. 5 and insist on fitting them 
with a power form in Eq. (|3TD , we would find that one can fit the last 12 mass points 
with 5 = 0.054 ± 0.028 and x /DF = 0.98. This is plotted as the solid line in Fig. 5. 
This serves as a caveat to show that one can misconstrue the contamination of the 
zero modes as the quenched chiral log. 

Finally, we look for the chiral log in fp, which according to Eq. (^) should grow 
in the form of (^)'^- We calculate the renormalized fp from 

/pa^ = lim ZpJGpp{t)2mpae"'^'/\ (34) 

t/a>l * 

where we insert the fitted mpa to obtain fpo^. The ^JCppit) is understood to have 
included the correction of the boundary reflection discussed in Sec. ^ in association 
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20^ Lattice with Overlap Fermions 




Figure 11: Renormalized fpo^/ f-^a vs m^a^ . The solid line is a fit over the lowest 14 
points. 

with the determination of f-j^. The renormalization constant Zp is the same as the 
scalar renormalization constant Zs- The latter is determined from the matching of 
renormalization group invariant quark masses at fixed pseudoscalar mass |]^ and the 
details will be given elsewhere [^. We plot the renormalized fpo^ in Fig. 10 as a 



function of moct. We see that it rises sharply at small quark mass. Since the log form 
in Eq. (|32| ) falls off slower at larger quark mass, it gives a more stable fit covering a 
larger range of the quark mass than the exponential form. Based on this observation, 
we fit fp in the log form 

fpa^ = fpa^{l - S[\n{Amoa/Ala^) + 1]} + Bmoa, (35) 



where we input A from the fit to mp in Eq. (|3^). It turns out that fits covering the 
whole range of 16 points have /DF greater than unity mainly due to the fact the 
fpo? of the two smallest quark masses are too high. Since errors of these two points 
are large and they are more susceptible to systematic errors from the boundary effect 
for such low pion masses, we shall ignore them and fit Eq. P5| ) with the other 14 
points. Again, results with reasonably small errors and x^/DF less than unity are 
reported in Table 6. 

We see that 5 obtained from fpo^ is consistent with those determined from m^a^. 
Combining all the fits from rn^o? and /pa^, 5 appears to span the range from 0.15 
to 0.4 which is consistent with that predicted from the quenched rj' loop in the chiral 
perturbation theory in Eq. (|27|) ||35|, ESl, although it tends to be on the high side. We 
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Table 6: fpo^, 6 and 'x^ /DF as fitted for fpo^ in Eq. ( |35D are listed. 



Ax (GeV) 


fpa' 


6 


B 


XVDF 


0.6 


0.083(2) 


0.35(4) 


0.36(2) 


0.93 


0.8 


0.060(3) 


0.48(7) 


0.36(2) 


0.93 



note that the general behavior of our data is close to what are observed in Ref. j^S], ^ 
but our 5 is about three to six times larger than theirs. 

To alleviate the apprehension that the apparent singular behavior of fpo^ may be 
caused by certain unknown boundary effect on the source in our lattice with fixed 
boundary condition, we consider the ratio of fp in Eq. (|3^ and /^r in Eq. (pOf ) 



r 2ir ,. ZpV2mpaGQ^A^p{p = Q,t) 
fpa Ua = hm — — — . 36 

t-—*oc 2moa Ga^p{p = 0, t) 

This should cancel out the boundary effect on the source. We plot the ratio fpa^/ fn^ 
as a function of m^a^. We see from Fig. 11 that the singular behavior is still visible 
in fpo?/ f^^a which confirms that the chiral log singularity is indeed present in fp. 



6 Summary 

To conclude, we have studied the chiral properties of the pseudoscalar meson on a 
quenched lattice with overlap fermions. The lattice size is 20^ with lattice spacing 
a = 0.148 fm set by the pion decay constant fj^. This gives a physical size of 3.0 fm 
which is about 4 times the Compton wavelength of the lowest mass pion. 

We first clarified the role of the zero modes in the pseudoscalar meson propagator 
in association with the generalized Gell-Mann-Oakes-Renner relation. We find that 
the zero mode contribution to the pseudoscalar meson propagators extends to a fairly 
long distance in the time separation and it is imperative to avoid it since it is a finite 
volume effect. Otherwise, one might be led to wrong conclusions about the behavior 
of the pion mass and decay constants near the chiral limit. For example, if one fits 
the pion mass by choosing a time window which straddles over the kink in Fig. 2, 
one will find that the pion mass squared does not approach zero at the chiral limit 
with a linear extrapolation in the quark mass. Interpreting it as due to the quenched 
chiral log, one may fit it with a chiral log form and obtain a positive 6 which can be 
misleading. In view of the fact that the quenched chiral log fit of m^a^ from data 
obtained from the shorter time range of Gpp{p = 0, t) as plotted in Fig. 5 has a much 
smaller 6 than the present study and is consistent with that obtained in Ref. it 
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would be useful to clarify that the shifted real modes (the 'would be' zero modes) in 
the Wilson action are not responsible for the observed quenched chiral logs in m^a^ 
and /pa^ in this case. 

Due to the chiral symmetry of the overlap fermion, we obtain the non-perturbative 
renormalization constant = 1.589(4) from the axial Ward identity. We find it to 
be fairly independent of mga which is an indication that the O(a^) error is small. The 
renormalized pion decay constant f^^ has a positive slope with respect to mga. With 
a small error (1.6%) and devoid of the complication of the quenched chiral log, is 
an ideal physical observable to set the scale of the lattice. 

We studied extensively the issue of quenched chiral log. We have clearly seen 
the quenched chiral log both in the pseudoscalar meson masses with equal quark 
masses and their pseudoscalar decay constants. Various fits give S to be in the range 
of 0.15 to 0.4 which is in accord with that estimated from the double r) propagator 
approximation of the quenched rj' loop which predicts it to be ~ 0.18. 

We finally have a reliable tool in the overlap fermion to study the chiral symmetry 
properties of hadrons at low energies including the quenched chiral logs. One should 
go to different lattice spacings to study the continuum limit in the future. 

This work is partially supported by DOE Grants DE-FG05-84ER40154 and DE- 
FG02-95ER40907. We wish to acknowledge enlightening discussions and exchanges 
with L. Giusti, M. Golterman, P. Hernandez, K. Jansen, D. Lin, M. Liischer, R. 
Sommer, A. Soni, and C. Rebbi. 
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